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ABSTRACT 

High cadence, multiwavelength, optical observations of solar magnetic bright points, captured at 
disk center using the ROSA and IBIS imaging systems on the Dunn Solar Telescope, are presented. 
Magnetic bright points manifesting in the Na i Di core are found to preferentially exist in regions 
containing strong downflows, in addition to co-spatial underlying photospheric magnetic field concen- 
trations. Downdrafts within Na i Di bright points exhibit speeds of up to 7 kms~^, with preferred 
structural symmetry in intensity, magnetic field and velocity profiles about the bright point center. 
Excess intensities associated with G-band and Ca ii K observations of magnetic bright points reveal a 
power-law trend when plotted as a function of magnetic flux density. However, Na i Di observations 
of the same magnetic features indicate an intensity plateau at weak magnetic field strengths below 
«150 G, suggesting the presence of a two-component heating process; one which is primarily acoustic, 
the other predominantly magnetic. We suggest that this finding is related to the physical expansion 
of magnetic flux tubes, with weak field strengths (ss50 G) expanding by ^^76%, compared to a ^-^44% 
expansion when higher field strengths (sal50 G) are present. These observations provide the first 
experimental evidence of rapid downdrafts in Na i Di magnetic bright points, and reveal the nature 
of a previously unresolved intensity plateau associated with these structures. 

Subject headings: Sun: atmosphere — Sun: chromosphere — Sun: photosphere — Sun: surface 
magnetism 



1. INTRODUCTION 

In solar physics, the Na i Di absorption profile 
has long been a difficult line to study, and indeed 
model. Line profile fluctuations often result from a 
combination of Doppler wavelength shifts, p-mode os- 
cillations and other sm all-scale dynamical phenomena 
(jEdmonds fc Hsulll983D . Furthermore, temperature, mi- 
croturbulence and magnetic fiel d strength can affect the 
associated Na i Di line d epths (jAthav fc Canfieidlll969l : 
[slaughter fc WilsonlllQTl . The formation height of the 
Na I Di core has often b een th e most difficult to interpret. 
lAimanova fc GulvaevI (|1976f ) used slitless spectrograms 
to estimate a core formation height of I300-I700 km. 
This is abo ve the formation height of the Ca ii K core 
(«1200 km; iBeebe fc Johnson! 1 1 9691) . and overlaps wit h 
that of the Ha core («I500 km: IVernazza et"aLlll98fl ). 
However, observations indicate that Na i Di core im- 
ages have littl e in comrn o n wit h their Ca ii K or Ha 
counterparts. lEibe et all (|200I[ ) have utilized theoreti- 
cal Na I Di response functions to suggest that the for- 
mation height varies depending on the local tempera- 
ture and/or velocity. The authors estimate a Na i Di 
core formation height of <1000 km, based upon the 
non-local t hermodynamic e quilibrium radiative trans- 
fer code of iCarlssonI (119861 ) . This is cons i stent with 
the work of Hb mita' C1960D, iFinsterle et al.l (|2004D and 
iSimoniello et al.. (, 2008.1 . who suggest Na I Di line-core 

|d.jess@qub.ac.uk| 



emission originates from the upper photospheric layer 
(<800 km). 

Recent Na i Di observations and simulations have 
revealed that Na i Di core brightenings sample mag- 
netic field concentrat ions in the quiet solar network 
(jLeenaarts et a l."2010') . Through a comparison of simul- 
taneous MDI magnetograms with ground-based Na i D2 
observations, iCauzzi et al.l (|2000( ) were able to suggest 
that sodium bright points are coincident in position, 
size and shape with corresponding magnetic patches. 
However, the relatively poor spatial resolution of MDI 
(«1.2") limited the study to network bright points in 
excess of w3" photospheric diameter. Thus, it is im- 
perative to study these types of bright features with an 
imaging system capable of high spatial resolution. In this 
paper, we utilize a high-cadence multiwavelength data 
set to search for correlations linking small-scale Na i Di 
intensity enhancements to the underlying photospheric 
magnetic field. 

2. OBSERVATIONS 

The data presented here are part of an observing se- 
quence obtained during 13:46 ~ 14:47 UT on 2009 May 
28, with the Dunn Solar Telescope at Sacramento Peak, 
New Mexico. The Rapid Oscillat ions in the Solar Atmo- 
sphere (ROSA: rJess et al]|2010b[ ) six-camera system was 
employed in conjunction with the Interferome tric Bldi- 
mensional Spectrometer fIBIS: ' Cavallini|[2006D to image 
a 69.3" X 69.1" region positioned at solar disk center. 
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TABLE 1 

ROSA/IBIS FILTER AND CADENCE OVERVIEW. 



Filter Used 


Exposure Time 


Frames per 


Reconstructed 




(ms) 


Second 


Cadence (s) 


G-band 


15 


30.3 


0.528 


Blue continuum (4170 A) 


10 


30.3 


0.528 


Ca II K core 


200 


3.8 


4.224 


Ha core (Zeiss) 


240 


3.8 


4.224 


LCP"^ (6302.375 A) 


240 


3.8 


4.224 


RCP'' (6302.375 A) 


240 


3.8 


4.224 


Na I Di (5895.94 A) 


70 


2.3^= 


39.703'* 



": Left-hand Circularly Polarized light 
*"; Right-hand Circularly Polarized light 
Average frames per second including Fabry-Perot tuning time 
"*: Cadence of a complete Na I Di profile scan 



IBIS sampled the Na i Di absorption line at 5895.94 A, 
incorporating nine wavelength steps with ten exposures 
per step to assist image reconstruction. More details of 
our instrumentation setup are given in table [T] A spa- 
tial sampling of 0.069" per pixel was used for the ROSA 
cameras, to match the telescope's diffraction-limited res- 
olution in the blue continuum to that of the CCD. This 
results in images obtained at longer wavelengths being 
slightly oversampled. However, this was deemed desir- 
able to keep the dimensions of the field of view the same 
for all ROSA cameras. IBIS employed a spatial sampling 
of 0.083" per pixel, which allowed ROSA's near square 
field-of-vicw to be contained within the circular aperture 
provided by IBIS. 

During the o bservations, high-order adaptive optics 
(jRimmeld I2004D were used to correct wavefront defor- 
mations in real-time. The acquired images were fur- 
ther improved through speckle reconstruction algorithms 
dWeigeh & Wirnitzcr 1983; Wogcr ct al. 2008), utihzing 
16 ^^ 1 and 10 -^ 1 restorations for ROSA and IBIS 
images, respectively. Post-reconstruction cadences for 
ROSA and IBIS are displayed in the fourth column of 
table [TJ A full image-reconstructed IBIS scan through 
the Na i Di absorption line resulted in a cadence of 
39.703 s, and includes a blueshift correction required due 
to th e use of classical etalon mountings (jCauzzi et alj 
|2008() . To insure accurate coalignment in all band- 
passes, broadband time series were Fourier co-registered 
and de-stretched using a 40 x 40 grid, equ ating to a 
» 1. 7" separation between spatial samples Oe ss et all 
120071 12008). Narrowband images, including those from 
IBIS, were corrected by applying destretching vectors 
estab lished from simultane ous broadband reference im- 
ages (jReardon fc Cavallini| [2008). Line-of-sight magne- 
tograms were generated as a difference image, normal- 
ized to their sum, of left- and right-hand circularly po- 
larized light, obtained 125 mA into the blue wing of the 
magnetically-sensitive Fe i absorption line at 6302.5 A. 
Utilizing high-resolution MDI magnetograms of a small 
active region, obtained in conjunction with ROSA us- 
ing an i dentical optical s etup, and following the previous 
work of llshikawa et all (j2007f) . these normalized differ- 
ence images were subsequently scaled to provide cali- 
brated line-of-sight magnetograms. Due to the ROSA 
cameras remaining within their linear regime, and MDI 
observations residing below their saturation threshold. 



only linear corrections were requi red to produce c ali- 
brated line-of-sight magnetograms ()Jess et al.ll2010cD . 

3. ANALYSIS AND DISCUSSION 

No active regions were present on disk during our ob- 
servations, and the ROSA/IBIS field-of-view contained 
no large-scale magnetic activity. However, a small col- 
lection of magnetic bright points (MBPs) were visible at 
heliocentric co-ordinates (3.9",— 1.5"), providing an op- 
portunity to examine these kiloGauss structures without 
the line-of-sight effects associated with near-limb obser- 
vations. Figure [T] shows snapshots of all the ROSA/IBIS 
sub-images obtained at 14:06:21 UT, and reveals a collec- 
tion of intensity enhancements, visible in the core of the 
Na I Di line (lower- middle panel of Fig. [T]). Ca ii K obser- 
vations were acquired through a 1 A wide filter, with re- 
sulting emission originating from the upper photosphere 
to the lower chromosphere (see e.g. lower-right panel of 
Fig- HI)- The Na i Di core image is an intensity map, cre- 
ated by establishing the line-profile minimum per pixel. 
By displaying Doppler-compensated line-center intensi- 
ties, rather than rest-wavelength line-center intensities, 
brightness variations throughout the image are more in- 
dicative of the source function than of the velo cities 
present in the line-forming region (jLeenaarts et al.12 010'). 
Examination of G-band (upper-left panel of Fig. [T]) and 
Na I Di core images reveal, crucially, that not all Na i Di 
core brightenings can be directly linked to underlying 
MBPs. Thus, there may be two possible interpreta- 
tions which promote efficient Na i Di core brighten- 
ing. The first implies independent heating processes, 
whereby magnetic and non-magnetic brightenings are 
produced via unrelated mechanisms. Contrarily, another 
interpretation revolves around the same heating mecha- 
nism, but that the process is somehow more effective in 
highly magnetic areas, such as those found in MBP lo- 
cations. Therefore, it is important to investigate small- 
scale MBPs, and establish which process allows enhance- 
ment of the Na i Di line above these < 2" G-band struc- 
tures. 

Information contained within the Na i Di core veloc- 
ity maps, generated through Doppler-shifts of the profile 
minimum, begins to unlock this mystery. An inspec- 
tion of Figure [1] reveals a preference for Na i Di and G- 
band MBP structures to exist in downflow regions (up 
to 7 kms~^). Furthermore, Ca ii K core and Ha core 
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Fig. 1. — Images from ROSA/IBIS obtained simultaneously at 14:06:21 UT. Color scales for the line-of-sight magnetogram and velocity 
maps indicate magnetic field strength (in Gauss) and Na I Di core velocity (in kms~^), respectively. Artificial downdraft saturation is 
displayed in the velocity map to aid identification of MBP features. Axes are solar heliocentric coordinates in arcseconds. Properties related 
to the MBP labelled 'A' are subsequently displayed in Figures [2] & [3] 



intensity enhancements are often observed to co-exist in 
locations of strong downflows (>3 kms^^) and longitu- 
dinal magnetic fields. Figure [2] shows a slice through a 
typical high red shift feature, displaying its G-band in- 
tensity, Na I Dl core intensity, magnetic field strength, 
and Na i Di line-core velocity as a function of distance 
across the underlying MBP. 

In the top-right panel of Figure [U we present an av- 
eraged Na I Dl line profile (solid line) from the central 
region of a MBP, labelled 'A' in FigurelU where high red- 
shifts are observed. A dashed line represents a Na i Di 
profile, which has been spatially, and temporally, aver- 
aged over the entire field-of-view and time series (in ex- 
cess of 9 X 10^ individual profiles). It is evident that 
the MBP profile has been significantly Doppler shifted 
when compared to the average Na i Di profile. The 
difference between these two profiles is displayed in the 
bottom-right panel of Figure [H and reveals a blue bump 
at w— 4 kms^^, alongside a red shifted dip at w8 kms~^. 
These asymmetries are indicative of shock formation 
below the chromospheric canopy, and are in remark- 
able a greement with the simulations of iLeenaarts et aU 



iRoudier et al.l ()2006[ ) investigated velocity flows asso- 
ciated with the Na i Di line in strong magnetic field 
regions surrounding active regions and sunspots. These 
authors detected an abundance of downfiows in these lo- 
cations, with velocities nearing 2 kms^^. As shown in 
the lower-right panel of Figure [21 the Na i Di core ve- 
locities presented here are considerably higher, reaching 
«7 kms~^ at the center of the MBP. Even artificially bin- 
ning o ur data to mimic the resolution of IRoudier et al.l 
(|2006D would not reduce the peak Na i Di Doppler ve- 



locities to below 3 kms ^. However, the spectral res- 
olution of IBIS (23 mA FWHM) is considerably better 
than the resolution obtained by the Multichannel Sub- 
tra ctive Double P a ss sys tem (144 mA FWHM) employed 
by IRoudier et al.l (|2006f ) , and allows a higher degree of 
precision when determining line-core velocities. Quanti- 
tatively, the Dop pler velocity error, err, c an be defined as 
err = ± fwhm c (|Langangen et al.ll2008l) . where A is the 
line-core wavelength and c is the speed of light. This pro- 
duces an associated IBIS velocity error of ±0.5 kms~^, 
compared to a n erro r of ±3.7 kms~^ in the results of 
IRoudier et"an (|2006f) . In addition, Na i Di downdraft 
velocities approaching 7 kms~^ are consistent wi th the 
magne to- hydrodynamic simulations of Lcenaart s et al.l 
(j201Clf ). Thus, for the first time, the combined high 
spatial and spectral resolutions of speckle-reconstructed 
Na I Dl IBIS data sets allow fine-scale Doppler velocities 
to be evaluated with an unprecedented degree of preci- 
sion. 

From the lower panels of Figure [21 we see that MBP 
structures found in Na i Di core imaging demonstrate an 
increased width over simultaneous G-band observations. 
This is consisten t with the "aureoles" , or halos, found by 
ILeenaarts et al.l ()201C1[ ). who interpret these additional 
brightenings as an artifact of strong radiative scattering. 
It must be noted that the decrease in normalized inten- 
sity to values below 1.0, at the edges of the G-band MBP, 
is due to intensity sampling of the dark inter-granular 
lane, where the MBP resides. Examination of Figure [H 
suggests that the fastest downdrafts are co-spatial with 
the center of the MBP, and that the velocities drop to 
approximately zero at the edge of the structure. Further- 
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Fig. 2. — The left-hand panels show the G-band (upper) and Na I Di core (lower) intensity images of the MBP labelled 'A' in Figure [Tl 
Overplotted in red are the pixels used to determine the excess intensity plotted in Figure|4] The remaining panels plot the one-dimensional 
variations in normalized G-band intensity (upper-middle), normalized Na I Di core intensity (lower-middle), magnetic field strength (upper- 
right), and Na I Di core downdraft velocities (lower-right), plotted as a function of distance from the center of the MBP along the cut 
indicated by the horizontal dashed lines in the left-hand panels. Green shading in the normalized Na I Di core intensity plot indicate the 
extent of the halo around the MBP, while red shading marks the locations of the central MBP. Blue shading reveals regions of non-zero 
downdraft velocities, synonymous with Na I Di core MBP enhancements. Cross-hatched regions display the values used when determining 
the excess intensity. 



more, the strongest longitudinal magnetic field strength 
is also co-spatial with the center of the MBP, and corre- 
lates well in size with the width of the G-band structure. 
Following analysis of all MBPs observed with red shifts 
exceeding 3 kms~^ and underlying magnetic field con- 
centrations, prevailing symmetry is found in intensity, 
velocity and magnetic-field strength, as a function of ra- 
dial distance from the center of the MBP. This indicates 
the preference for an ordered structure, whereby maxi- 
mum values are attributed to the center of the MBP, and 
steeply decaying measurements are found at increasing 
distances from the center. 

Investigating the intensity excess of MBPs, which cor- 
respond to all MBPs containing strong (>3 kms~^) 
downflows and G-band inter-granular brightenings, re- 
veals a correlation with the underlying magnetic field 
concentrations. Figure U reveals how the intensity ex- 
cess of G-band, Ca ii K and Na i Di observations varies 
as a function of the magnetic flux density, $, defined as 
the total magnetic field strength divided by the num- 
ber of occupied pixels. The excess intensity. A/, is de- 
termined as an averaged pixel intensity for the entire 
MBP, following a subtraction of intensities equal to the 
modal value observed in quiet areas (central panels of 
Figure [5]). For G-band images, MBPs are isolated from 
granules a nd inter-granular lane s using the automated 
routines of iCrockett et al.l ()2009f ). These isolated pixels 
also provide the structures used to determine the mag- 
netic flux density from the line-of-sight magnetograms. 
Na I Di core MBP observations consist of a bright central 
structure, surrounde d by a dimmer halo, r esulting from 
radiative scattering (jLeenaarts et al.ll2010[ ). The central 
MBP is isolated from the halo by utilizing simultane- 
ous Na I Di core velocity maps. The lower-right panel of 
Figure [2] reveals how the downdraft velocities drop to ap- 



proximately zero at the edge of the bright central MBP. 
Therefore, pixels where the downdraft velocity is non- 
zero are used to isolate the central MBP regions, thus 
minimizing contributions to the excess intensity caused 
by radiative scattering, which is only weakly dependent 
on the magnetic field strength. 

In the case of G-band observations (upper panel of 
Fig. m, a linear depende nce observe d in t he log-log 
plot is in agreement with iBeck et all (|2007[ ). who de- 
tect identical trends for G-band brightenings in ac- 
tive region moats. Through analyses of network bright 
points (>3" in diameter) and large-scale active r egion 
structures (>1 5 " in di ameter), iGauzzi et al.l ()2000[) and 
iSchriiver et al.l ()1989aD were able to demonstrate a simi- 
lar trend for these features when observed in the Na I D2 
and Ca 11 K lines, respectively. Indeed, for the mag- 
netic fiux densities presented here (above wl50 G in the 
case of sodium), the calcium, G-band and sodium ex- 
cesses (Fig. lU are best fitted by a power law of the type 
A/ oc $^, where /3 is the power-law exponent. Here, we 
determine (3 = 0.61 ± 0.02 for Ca 11 K observations, and 
(3 — 0.59 ±0.03 in the case of Na i Di data. These values 
of power-law exponent are in agree ment with the pre- 
vious Na I D2 study of IGauzzi et al.l (|2000. ) . and remain 
consistent with exponeii t s found in Ga 11 K inve stigations 
(ISkumanich et all [TOTSt fNindos fc ZirmI llQQSl) . While 
these consistencies with previously published work may 
seem initially superfluous, the spatial scale over which 
these results agree is extraordinary. Schrijver et al.l 
(|1989aj) investigated large-scale magnetic structures con- 
tained within an active region, often with spatial sizes 
approaching 70" in diameter. The MBPs investigated 
here consist of structures as small as ?a0.2" in diameter, 
more than two orders of magnitude smaller than those 
studied previously. 
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Fig. 3. — Top left: G-band intensity image of the MBP labelled 'A' in Figure[T] Lower left: Line-of-sight velocity map of the same MBP, 
obtained from Doppler wavelength shifts of the Na I Di profile. The color scale indicates flow velocities in kms~^, where red and blue 
colors represent downflows and upflows, respectively. In each panel, the solid black contour indicates regions where downflow velocities 
exceed 4 kms~^. Top right: Averaged Na I Di profile (solid line) generated from locations contained within the contours displayed in the 
left panels of this Figure, overplotted with a rest-wavelength Na I Di proflle (dashed line). Lower right: The difference between the two 
proflles. In both right-hand panels, the solid vertical line represents the Na I Di rest wavelength at 5895.94 A. 



At magnetic flux densities below sal 50 G, there is an 
appreciable "turn-off" from the power-law fit applied 
to the Na i Di observations. The lower panel of Fig- 
ure |4] reveals how the Na i Di intensity remains elevated, 
even as the magnetic flux density begins t o diminish- 
Throug h analysis of Ca ii K observations. Schriiver et al.l 
(|1989a| ) suggest large deviations in line-core intensity, at 
relatively small magnetic flux densities, may be a re- 
sult of a large inclination of the magnetic flux tubes, 
or small-scale energetic phenomena. Inclination effects 
on the data presented here are not believed to be sig- 
niflcant for two reasons. Firstly, the inclination an- 
gle of MBPs, with respect to the vertical, is typicall y 
less than 10° ([Sanchez Almeida fc Martinez Pilleg|l994D . 
Secondly, due to the observations being acquired at disk 
center, line-of-sight effects are minimized. This allows 
photospheric structures to remain co-spatial when ob- 
served in chromospheric passbands, and since no feature 
migration is found between G-band and Ca ii K atmo- 
spheric heights, further solidifies the presence of minimal 
inclination angles. Additionally, the intensities of MBPs 
corresponding to G-band, Ca ii K and Na i Di observa- 
tions show no temporal v ariations that may be suggestive 
of microflare ([Jess et al.l l 2010aD and/or EUerman bomb 
(jEllermanlll917[) activity. Through examination of pho- 
tospheric G-band MBPs, we observe the physical area 
of such structures to reduce as the magnetic fleld be- 
comes weaker, corrob orating the observational results of 
llshikawa et al.l (2007). Contrarily, we find MBPs, which 
correspond to magnetic flux densities below the turn- 
off at Ril50 G, show similar two-dimensional areas when 



viewed in the Na i Di line core. An average Na i Di core 
area of O.lliO.Ol Mm^ (47±4 pixels) is found for these 
structures, suggesting their excess intensity plateau be- 
low «150 G may be r elated to the physica l expansion of 
magnetic flux tubes. ^Sol anki et al.l (jl999f ) utilize mod- 
els based upon infrared spectral lines to infer how small 
magnetic flux tubes expand more rapidly as a function of 
atmospheric height, when compared to larger magnetic 
structures. Thus, in order to maintain relatively constant 
physical sizes present in Na i Di observations, it is nec- 
essary for the physically smaller (and weaker) magnetic 
fleld concentrations in the photosphere to expand at a 
greater rate. To quantify this result, we find MBPs at 
the lowest magnetic flux densities («50 G) show a ~76% 
increase in area between G-band and Na i Di observa- 
tions, compared to a ~44% increase for MBPs near the 
«150 G turn-off. 

A discontinuity in the correlation between the mag- 
netic field and Na i Di intensity, for magnetic flux den- 
sities below «150 G, suggests that these MBP inten- 
sities are independent of the underlying magnetic fleld 
strength. Combining N a i Di obs ervations with 3D non- 
LTE simulations, Lcen aarts et a l. (2010) conclude that 
the Na i Di line is not chromospheric in nature, and 
instead samples the atmospheric layer below the mag- 
netic canopy, where magneto-acoustic shocks are believed 
to dominate. Our findings may therefore indicate a 
two-component heating process; one which is primarily 
acoustic and dominates the heating of low magnetic field 
strength areas, while the other process is predominantly 
magnetic in nature. Within the presented data, MBPs 
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Magnetic Flux Density (Gauss) 

Fig. 4. — Intensity excess for each MBP containing downdraft ve- 
locities exceeding 3 kms"^, plotted as a function of MBP-averaged 
magnetic flux density for G-band (top), Ca II K (middle) and 
Na I Di (bottom) observations, and displayed on logarithmic axes. 
In each panel, power-law fits have been applied to the data which 
best minimizes the chi-squared error statistic. 

with velocities exceeding 3 kms"^ represent the dominant 
population. However, as MBPs can exist with velocities 
below this threshold, their excess intensity behaviour, as 
a function of magnetic flux density, may differ as a re- 
sult. We note that two-component heating models have 



previously been proposed to interpret observed chromo- 
spheric flux densities in the most inactive solar-type s tars 
(jSchriiver et al.lll989bHMathioudakis fc Dovlelll992[ ). 



4. CONCLUDING REMARKS 

Joint observations, acquired with the Rapid Oscilla- 
tions in the Solar Atmosphere (ROSA) and Interfero- 
metric BIdimensional Spectrometer (IBIS) imaging sys- 
tems, have enabled unprecedented views of solar mag- 
netic bright points (MBPs). These data, obtained si- 
multaneously with high spatial, temporal and spectral 
resolutions, have revealed how MBPs manifesting in the 
Na I Di core are found to preferentially exist in regions 
containing downflows, in addition to co-spatial under- 
lying magnetic field concentrations. Downdrafts found 
within Na i Di bright points exhibit speeds of up to 
7 kms^^, with preferred structural symmetry in intensity, 
magnetic field and velocity profiles about the center of 
the MBP. Excess intensities associated with G-band and 
Ca II K observations of MBPs demonstrating downdrafts 
exceeding 3 kms~^ reveal a power-law trend when plotted 
as a function of magnetic flux density. However, Na I Di 
observations of the same magnetic features indicate an 
intensity plateau at weak magnetic field strengths below 
«150 G. We interpret this phenomena as the physical ex- 
pansion of magnetic fiux tubes, with weak field strengths 
(«50 G) expanding by ^76%, compared to ~44% ex- 
pansion for structures at higher magnetic field strengths 
(«150 G). 
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